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Abstract 



< 

(— i | We consider a conforming finite element approximation of the Reissner-Mindlin 

system. We propose a new robust a posteriori error estimator based on H(div ) con- 
forming finite elements and equilibrated fluxes. It is shown that this estimator gives 
rise to an upper bound where the constant is one up to higher order terms. Lower 
bounds can also be established with constants depending on the shape regularity of 
the mesh. The reliability and efficiency of the proposed estimator are confirmed by 
some numerical tests. 

t> 
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AMS (MOS) subject classification 74K20, 65M60, 65M15, 65M50. 

1 Introduction 

^ ■ The finite element method is often used for the numerical approximation of partial differ- 
ential equations, see, e.g., [7J El EES]. In many engineering applications, adaptive techniques 
based on a posteriori error estimators have become an indispensable tool to obtain reliable 
results. Nowadays there exists a vast amount of literature on locally defined a posteri- 
ori error estimators for problems in structural mechanics. We refer to the monographs 
[H [21 [291 E2] for a good overview on this topic. In general, upper and lower bounds are es- 
tablished in order to guarantee the reliability and the efficiency of the proposed estimator. 
Most of the existing approaches involve constants depending on the shape regularity of the 
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elements; but these dependencies are often not given. Only a small number of approaches 
gives rise to estimates with explicit constants, see, e.g., [TJ El [151 EH IZEJ |25j [281 [211 130] . 
However in practical applications the knowledge of such constants is of great importance, 
especially for adapt ivity. 

The finite element approximation of the Reissner-Mindlin system recently became an 
active subject of research due to its practical importance and its non trivial challenges to 
overcome. In particular, appropriated finite elements have to be used in order to avoid shear 
locking. Such elements are in our days well known and different a priori error estimates 
are available in the literature. On the contrary for a posteriori error analysis only a small 
number of results exists, we refer to [51 El HH [121 El [261 [23 [21] . Most of these papers enter 
in the first category mentioned before and to our knowledge only the paper [21] proposes 
an estimator where an upper bound is proved with a constant 1. Hence our goal is to 
give an estimator that is robust with respect to the thickness parameter t, with an explicit 
constant in the upper bound, that is also efficient and that is explicitly computable. For 
these purposes we use an approach based on equilibrated fluxes and if(div )-conforming 
elements. Similar ideas can be found, e.g., in [61 [151 ETJ [281 EH]- For an overview on 
equilibration techniques, we refer to [H 125] . 

The outline of the paper is as follows: We recall, in Section 2, the Reissner-Mindlin 
system, its numerical approximation and introduce some useful quantities. Section 3 is 
devoted to some preliminary results in order to prove the upper bound. This one directly 
follows from these considerations and is given in full details in section 4. The lower bound 
developped in section [5] relies on suitable norm equivalences and by using appropriated 
H(div) approximations of the solutions. Finally some numerical tests are presented in 
section El that confirm the reliability and the efficiency of our error estimator. 



2 The boundary value problem and its discretization 

Let £1 be a bounded open domain of M 2 with a Lipschitz boundary T that we suppose to be 
polygonal. We consider the following Reissner-Mindlin problem : Given g G L 2 (Q) defined 
as the scaled transverse loading function and t a fixed positive real number that represents 
the thickness of the plate, find (oo,4>) G Hq(Q) x Hq(Q) 2 such that 

a(4>,i>) + ( 7 , Vv-xP) = (g,v) for all (v^) G H^Q) x H^Q) 2 , (1) 

where 

7 = At" 2 (Vu - 0) and a(<j>, if)) = / Ce((j))e(<ip)dx. (2) 

Jn 

Here, (• , ■) stands for the usual inner product in (any power of) L 2 (Q), the operator : 
denotes the usual term-by term tensor product and 

£(0) = ^(V0+(V0f). 
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C is the usual elasticity tensor given by 

Ce(<f>) = 2 n e(<p) + A tr(e(4>)) 1. 

The parameters /i, A and A are some Lame coefficients defined according to the Young 
modulus E and the Poisson coefficient v of the material. In the following, for shortness the 
L 2 (Z))-norm is denoted by || ■ \d. The usual norm and seminorm of H l (D) are respectively 
denoted by || • and | - |i,r> and the usual norm on is denoted || ■ ||-i,z> For all 

these norms, in the case D — Q, the index Q is dropped. The usual Poincare-Friedrichs 
constant in Q is the smallest positive constant cp such that 

||0||<c F |0|i V0G# o W- 

By Korn's inequality [22], a is an inner product on Hq(Q) 2 equivalent to the usual one. 
Indeed, defining the energy norm 1 1 ■ | \c by 

it can be shown (see annex IT. 1 h that 

H?<-|HI^G# W- (3) 
A 4 

Consequently, the continuous problem ([I])-© is well-posed. 

Lemma 2.1 The problem (QP-^ has a unique solution (u, 0) G Hq(Q) x Hq(Q) 2 . 

Proof: Defining the functional F((u, 0), (v, if))) = a(0, i/))+('j, Vv—tp) with 7 = A t~ 2 (Vw- 
0), let us establish its coerciveness, namely that there exists k > such that 

F((cu, 0), [cu, 0]) > k (\co\l + |0| 2 ), V(w, 0) G H*{n) x H^Q) 2 . (4) 

Fix an arbitrary pair (oj, 0) G Hq x Hq(Q) 2 . First of all, and the standard Cauchy- 
Schwarz inequality lead to 

F((u;,0),(u;,0)) > /i|0| 2 + At" 2 ( (1 - V )\cu\l + (l - -) ||0|| 2 V Vr? > 0. 
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Then we directly obtain 



F((u, 0), ( Wj 0)) > ^|0| 2 + Ar 2 (l - r^l 2 + + xr 2 (l - I ) ) ||0|| 2 , V77 > 0. (5) 



2c|At 
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Choosing now 77 = ^— — = < 1 in ©, we have 

/i + 2c z F Xt~ 2 



i- 



This shows that flU) holds with k = min ( — , — - „ - — - ). The conclusion follows from 

the Lax-Milgram lemma for which the other assumptions to fulfill are obvious. ■ 



Let us now consider a discretization of (DO)-© based on a conforming triangulation Th 
of Q composed of triangles. We assume that this triangulation is regular, i.e., for any 
element T G Th, the ratio h?/ Pt is bounded by a constant a > independent of T and 

of the mesh size h = max/i^, where is the diameter of T and the diameter of its 

TeTh 

largest inscribed ball. We consider on this triangulation the classical conforming Pi finite 
element spaces Wh x ®h defined by 

W h = {v h G C (ti);v h = 0ondn and v h \ T G ¥ 1 (T) VT G T h } C H^O), 

B h = W h xW h C H^Q) x H^Q). 

The discrete formulation of the Reissner-Mindlin problem is now to find (uh, (j>h) £ Wh x 0^ 
such that 

iph) + (7fc, - = 0, v h ) for all (u h , <E W h x Q h , (6) 

with 

7 fc = Ar 2 (Vw h -IW h ). (7) 

Here, R^ denotes the reduction integration operator in the context of shear-locking with 
values in the so-called discrete shear force space Th which depends on the finite element 
involved [31 HJ [18j [19j [31] . We assume moreover that 

For all tj) h G Q h ,R h ^ h G H (rot,Q), 

where H (rot, Q) = {v G L 2 {Vt) 2 ] rot v G L 2 {Vt) and f • r = on <9f2}, equipped with the 
norm 

IMIjr(rot,n) = IMIo+ ll ro Hln- 

Here, for any v = {vi,v 2 ) T G L 2 (Q) 2 , rotv = dv2/dx — dvi/dy and r is the unit tangent 
vector along dQ. In this work, R^ is defined as the interpolation operator from 0^ on the 
Ho(rot,Q) conforming lower-order Nedelec finite element space [22] . 

By the usual Helmholtz decomposition of any H (rot,Q) vector field [8j p. 299], there 
exists w G H^(Q) and G H^(Q) 2 such that : 

(R„ - I)<f> h = Vw-P, (8) 

as well as a constant C > such that 

H{rot,n)- 
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More precisely, we introduce the constant cr such that 



W < c R \\rot(R h - I)<f> h \\, 

which can be evaluated by [22] 

( . , (divv,q)^ 

Cr = ml SUp 

\^LHn) veH i (n)2 \\q\\ |u|i 

Given the exact solution (w,0) G Hq(Q) x ^(fi) 2 as well as the approximated one 
(u>h, 4>h) G W/, x the usual error e£°* is defined as 

«) 2 = |w - co h \l + |0 - 4> h \l + A^ 2 || 7 - 7.f + A^V||rot( 7 - 7/Of + \\l ~ Tfcll-i- (9) 
The residuals are also defined as follows 

R eSl (v) = (g,v) - (7 hj Vu) for all v G ^(fi), (10) 
# eS2 (^) = -a{M) + ( 7A ,V) for all ^ G ^(fi) 2 . (11) 

Finally, let us now introduce, in the spirit of [21], the spaces Ndi V (Q) and Hdi V (fl) respec- 
tively defined by 

iWfi) = {y£ L 2 (n,R 2 )\dtvy G L 2 (ft)}, 
^dto(fi) = {xe L 2 (n,M 2 s )\divx G L 2 (fi,M 2 )}, 

where M. 2 S is the space of symmetric tensors of second rank. We now fix an arbitrary 
y* G Hdi v (Q) such that divy* = —Yl^g, where Uh is the projection operator from L 2 (Q) to 
the piecewise constant fonctions on the triangulation. Let us also fix x* G Ndi v (£i) such 
that div x* = —jh- Their existence and construction will be explained later on. 

We finally need to introduce the following mesh-dependent norm. For all (ip, v) G 
Hq(Q) x H^(n) 2 , we define 

lll(<MIII 2 ,,= IIW|| 2 + J2 whr\^ v ~ ^ ^ 12 ) 

reT h T 

For all functional F defined on Hq(TL) x Hq(Q) 2 , the dual norm associated with f fT2|) is 
classically defined by 

III Fill F ^ v ) hi\ 

\\W\\\-\,h= SUP — rrrr—. (13) 

(^)etf<J(f2)xtfi(fi) 2 \{0} III W v )\\\l,h 
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3 Preliminary results 



The aim of this section is to prove four lemmas which will be used in the following of the 
paper. 

Lemma 3.1 Let us consider (a,e) G (M^) 2 . Then we have 



X(r 2 -a 2 )\\V(u-cu h ) 



R h 4> h )\\ 2 + Xa 2 (l-2e)\\V(u-uj h )\\' 



- Rh<p, 



< A-V|| 7 - 7.II 2 - A« 2 (l - fj U h - R,M 2 -\a 2 (l-- e -e 

Proof: We first write 

||V(w - u) h ) - (0 - (f> h ) - (cf> h - RM\\ 2 

= - oo h )\\ 2 + ||0 - <j) h \\ 2 + Uh - R*0fc)|| 2 

-2(V(w -u h ),(j>- (p h ) - 2(V(w - u h ), 4> h - R h( f) h ) + 2(0 
Consequently, we have 

^Wl-lhf = A(r 2 -a 2 )||V(a;-a;,)-(0-R,0,)|| 2 

+ Aa 2 (|| V(w - w & )|| 2 + ||0 - 0,|| 2 + ||0, - R,0,|| 2 ) 
+ 2 Aa 2 (0 - ft , <p h - R h <p h ) - 2Aa 2 (V(w -u h ),<f>- <f> h ) 
- 2Aa 2 (V(w -u h ),(j> h - R h <Ph)- 
Using the three following Young inequalities 



-2(0 - 0^,0^ - R, 



< e\\<f>-<f> h \\ 2 + ^\\<i> h -Rh<j> h \ 

1„ 



2{V{u-u h ),(j)-(j) h ) < e\\V(u-cu h )\\ 2 + 



-||<^ ~ <Ph\\ , 
1„ . _ . „ 2 



2(V(w-u; h ),0 /i -R h (p h ) < s\\V(co-u h )\\ 2 + -\\(f) h -R h ^ h 

e 



we get 



A(r 2 - a 2 ) || V(u; - - (0 - R,0,)|| 2 

< A-4 2 || 7 - 7hf - Aa 2 (|| V(w - c^)|| 2 + ||0 - ft || 2 + ||0, - R,0,|| 2 ) 
+Aa 2 ^||0 - ft || 2 + -||0 h - R h 4> h f + e||V(w - c^)|| 2 + -||0 - 0J 2 



+ £ ||V(a;-a; h )|| 2 + -||0 ft -R h ,0 h 115 



\ „.2 



A"V|| 7 - lh f - \a z [I - ?J \\<j> h - R h <f) h \\ 2 -\a 2 (l-^-e 
-Xa 2 (l-2e)\\V(u-u h )\\ 2 . 



This proves the lemma. 
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Lemma 3.2 we have 

||7-7fc||-i <^^ + X)H-<t> h \\ 2 c + 2\\Res 2 \t v (14) 
Proof: First, it can be shown that for any ip G (Hq(Q)) 2 , 

\mi<2^+~xmi, 

so that 

(l-lh,i>) = a(<j> - </> h ,i/>) + a(<j> h ,if>) - (7ft, V) 
= a(0 - 0ft, ^) - Res 2 {ip) 

< ll^-^llcll^llc + ll^ll^l^li 

< ((2 (/i + \)) l/2 U - 0,11c + \\Res 2 \u) 

Hence we get 

||7-7j-i < ((2(/, + A)) 1 / 2 ||0-^|| c + || J Re S2 ||_ 1 ) 2 
<4(/i + A) ||0-0ft|| 2 + 2|| J Res 2 || 2 _ 1 . 

■ 

Lemma 3.3 

110 - 4>h\\c + A _1 t 2 ||7 - 7,|| 2 = i?esi(w - u h + w) + Res 2 (4> -</>h + P)- a(4> ~ </>h,P), 
where w and (3 are given by the Helmholtz decomposition 

Proof: This result is similar to the one given in [IT] . First, ([1]) and ([8]) lead to 

(7 - 7ft, (Rft - I)<f> h ) = (7 - 7ft, Vtu - P) 

= (7, Vw) - (7, p) - ( 7 ft, - 13) 
= (g, w) - a(0, 13) - (7ft, Vw; - (3) 
= -a(0 - 0ft, 13) + (#, iw) - a(0ft, /3) - (7,, Vw - 0). 
A simple calculation shows that 

||0-0ft|| c + A- 1 t 2 ||7-7/ l ir 

= a(0 - 0,, - <() h ) + (7 - 7,, (Vw - Vu h ) - (0 - 0,)) + (7 - 7/i > ( R /i - -00/0 
= (#, u - u h ) - a(0ft, - 0ft) - (7,, V(u - Wft)) 

+ (7ft, - 0ft) - a(0 - 0ft, 0) + (#, u; ) - a(0ft, /3) - (7/,, Vw - 0) 
= Res 2 (<p - (fi h + f3) + (g,u - u h + w) - (j h , V(u - u h + w)) - a(0 - h , (3) 
= Res 2 ((f) - 0, + (3) + Resi(u - u h + w) - a(0 - 0ft, /3). 
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So we get 

||0-<Mc + A- 1 t 2 || 7 -7, i || 2 = Re Sl (Lo - uu h + w) + Res 2 ((J) - <f) h + (3) - a((f) - (f) hl (3). 

■ 

Lemma 3.4 

lu-^ + /3\\ 2 c + lu-M\ 2 c + ^A-V|| 7 - lh \\ 2 
+ \ H ^rfj\\^^-^ + w)-(<p-<p h + f3)\\ 2 T 

Ten T 

< Re Sl (u -u h + w) + Res 2 (<p -<p h + f3) + -\\/3\\ 2 c . 

Proof: The proof is once again similar to the one in [IT]. Because of (IE]), we first remark 
that 

7 - 7h = Ar 2 (Vw - Voj h - + <p h + Vw - 
so that we have for all T G Th 

- cu h + w) - (<j> - cj> h + (3)\\t < A~ 2 t 4 ||7 - lhf T - 

Then, 

\U-<Ph + P\\ 2 c + \U-MI + \^t 2 h - luf 
A 



+ \ Yl t 2 + h 2 ~ u; h + w)-(<f)-<j) h + y 
TeT h T 

< - & + + \u - Mil + \^t 2 h - ih\\ 2 + \^ l t 2 £>- 7,11 

Ten 

< \-H 2 1| 7 - lh f + -a((f> -4> h + (3,4>-(P h + f3) + -a((f> -4> h ,4>- <f> h ) 

= A-^117 -lh\\ 2 + \{U- Ul + 2a(0 - 0,, /?) + \W C ) + \\\<P ~ Ml 
= \\<j> - Ml + A-VH7 " 7.II 2 + ^1 W + a(0 - <f> h , (3). 



From lemma 13.31 we get 

iM-fa + PWl + lM- Ml + ^A-V|| 7 - 7.II 2 

< Re Sl (uj - uj h + w) + i?es 2 (0 - h + /?) - a(0 - /3) + -||/3||^ + a(0 - h , /3) 

1 ^ 
= itesi(w - w/, + w) + i2es 2 (0 - <\>h + /?) + ^Pllc- 
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4 Reliability of the estimator 

Theorem 4.1 Let us consider < e < 1/2, as well as two parameters v\ > and v 2 > 0. 

Moreover, let us define 

/3 - + s — 1 
A(e) = max - + c|-^- — + 4(/x + A); 1 



/x >(l-2e) 



A(l-2e) / ' 



Then, 



f h ot ? < A 1 \ \\Re Sl \ \\\ h + A 2 ||i?e S2 ||^ + A 3 \\<f> -<p h + 

+A 4 \\<f )h - R h (j) h \\ 2 H{rot,a) ~ J2 A s\\ V ( w -Uh + w) 



with 



Ten 

A, = ^Aie) 2 ; 
A 2 = v 2 A{ef + 2; 

A 3 = -(- + -) -A(e); 

2 



A 4 = max 



- 1 



1 - 2e 
XA(e) 1 



2 + 2A(e)(^ + A)4 



+ 



(15) 



Proof: First of all, by using lemma [341 and the fact that < e < 1/2, we get 



i+izi 

At ^ ° F /x(l - 2e) 



A(l-2e) 



A-V|| 7 -7>iir 



H^-s-ll^ " R^ll 2 + A-^ 4 |M( 7 - 7 ,)lr + 117 - 7,|| 2 -i- 



1 - 2e 

Then, because of lemma 13.21 as well as 



we obtain 



D rot\2 



A -2 t 4|| rot(7 _ )|| 2 < _ ||0 _ faW* + 2 ||rot(0 h - R^,) 

A 4 



< I - + °F J^2^ + 4 (/^ + A) J ||0 - M\l + 2\\rot{<f> h - KM 
a \ / 2 



1 



A(l - 2e) 
+2||i?es 2 || 2 _ 1 . 



1 1 -D 

rt/, 



A-V|| T - 7 ,||^ + 



1 - 2s 



) h - r\, h (p h 



By the definition of A(e) as well as lemma I3.4[ we get 

(e r h ot ) 2 < A(e) (2Re Sl (uj - u h + w) + 2Res 2 {<j) - <\> h + (3) + 

-M-fa + PWl-Y, ^p-||V(a;-^ + «;)-(0-^ + /3)||2) 



2 1 



TeT h ' T 



1 - 2s 



+ ( -h^T- Hh - RhM 2 + SII^H 2 ! + 2\\rot(<j) h - R b ^ 



We notice that 

Res^u - u h + w) < |||i2eai|||_i jh |||(V',£*;-a; h + u;)||| 1 , h VV' e #o(^) 2 > 

Introducing now the parameters v\ > and ^2 > and using two times Young's inequality 
lead to 

(er) 2 < ^(6)111^111^+ i|||(^ W -a; fc + «;)|||? ifc 

+V2 A\e)\\Res 2 \\ 2 _ l + ^\ct ) -ct )h + P\l 
-A{e)U-<p h + P\\l + A(e)\\p\\l 

+ f fPP) Uh - Rh<f>hf + nRes2\\\ + 2||rot(0 ft - R^)|| 2 



1 - 2s 



E(^|)iiv(— ) 



Finally, choosing ?/> = — <ph + 0, we get 



\\\(^u-u h + w)\\\l h = ||v(0-^, + /3)|| 2 + £ —Lj-Hvcw-Wfc+uO 



, 2 , 72 II v v /* / \r *r/i My llT> 

and so ffT5l) holds. 



Corollary 4.2 Lei us assume that t < ^3\c 2 F /fi, and let us define : 

1 1 

y 2A 

T/ien, 



£ = max — , — J> . (16) 



2 < 2C f- + ^(3 + 2V3) + 4(/i + X)) I HiJeaxI || 2 1 
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2C + ^(3 + 2V3) + 4(// + A)^ + 2^ ||i?e S2 ||^ 

+ max ^7 + 4V3 ; 2 + ^ + ^(3 + 2^) + 4(/i + A)^ 2(/x + A)c^ ||<^ - ^^||i M , n) . 

Proof: Assuming 1 — 2e > 0, the parameters z/i and i/ 2 arising in the values of A 3 and 
Aj in ( II 5p are first chosen such that A 3 < and A J > VT G 7^. Namely we take 
v\ = v 2 = 2(/A(e). Consequently, we obtain 

(e^) 2 < iilll^lll^ + ^ll^ll^ + ^ll^-R^lll^n), (17) 

with 

' A x = 2(A(e); 
A 2 = 2(A(s) + 2; 

/ 2 - 1 _ 
A 4 = max f^r- ; 2 + 2A(e)(fi + A)4 

Now, in order to provide a result as sharp as possible, it remains to choose appropriately 
the parameter e to make the coefficients A±, A 2 and A 4 arising in ( IT7|) as small as possible. 
Since we always have 1 < 3/fi + 4(/i + A), the assumption t < ^/3Ac|,/yU leads to 

1 

At this stage we remark that the two functions A(e) as well as reach their minimum 

value for the same value of the argument e, namely for e = 2 — y/3. So, by a simple 
calculation, corollary 14.21 holds. ■ 

Now, it remains to bound each of the two residuals. 

Lemma 4.3 Let N G N* be such that maxF(T) < iV ; with Y(T) = #{T' G T h \ V C uj t } 

T&Th 

and ut = {K G T\K fl T 7^ 0} is t/ie paic/i 0/ elements surrounding T (consequence of the 
mesh regularity) . Then there exists k 2 > only depending on the mesh regularity such that 

llli^HI^ < 2Nk 2 2 J2 (t 2 + h 2 T )\\j h - y*f T + osc 2 (g), (18) 

TeT h 

where osc(g) corresponds to an oscillating term. 

Proof: For any v G Hq(Q), let us consider Vh = Jv where J : Hq(Q) — > Wh is defined such 
that (see, for example [H], known as the Clement operator) 

3ki>0; VTg %, \\Vv h \\ T < ^11 Vv\l T . (19) 
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Moreover, it can be shown [TT] that there exists k>2 > and K3 > such that for all T G Th 
and for any ip G H^(Q) 2 , 



\\V(v-v h )\\ T < KadlVu 



I W T 



MIWII 



CUT / 5 



/^t 1 !^ - v h \\ T < k 3 (\\Vv - ^|| WT + MIWLt) • 
Then for all v G flg(fi), we get 

Res\{v) = Resi(v — Vh) 

= (g,v- v h ) - (7 h , V(u - 

= (9 + divy*, v-v h )- (7 A - y*, V(f - u h )) 

TeT h 

< h T\Jt 2 + h 2 T \\g + divy*\\ T x T \\v - v h \\ T 

TeT h V + "T 

l|V(«-t; fc )|h 



So, we can write 



Resx(v) < ^ h T ^Jt 2 + h^Wg + divy*\\ T 
Ten 



x- 



< [ Yl 4h T (t 2 + h 2 T )\\g + dwy*\\ 2 T 



h T \\V^\\ 



+ MIW>|U) 



\TeT h 



1/2 



4(t 2 + h 2 T )\\ lh 

TeT h 



y*\\ 2 T 



1/2 



2 £(^iiv«-ve+ 
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Now recalling that max7(T) < N we have 
Ten 



( 



s <ny: 



TeT h 



t 2 + h 2 T 



<i 



Ten 

V TeT h l + 



^ <#III(MIIIU- 



So we get 



Resxiv) < I k\ ^ h T( t2 + ^r)ll# + divy 



* l|2 
T 



Consequently 



\ 1/2 

+ K 2 (t 2 + 4) || 7 * - y* III x v^iv| || (V, v) I Hi,, 

Ten J 
l-Reailll-i,/, <2iV (/eg h%{t 2 + h 2 T ) \\g + divy*f T 



Ten 
,*ii2 ^ r< a2 



Since dii;?/* = — IL^, we get ||g + divy*\\^ < C h\ \\g\\ 2 and ffTHl) holds. 

Lemma 4.4 For \I> G if^fi) 2 , we have 

Res 2 W < WC-^ix'-Ceifam U\\ e . 
Proof: Using standard Green formula, we easily obtain 

Res2(ip) = / (x* — Ce{4>h)) '■ s(i/j) dx + / (7^ + div x*) ip dx, 



(20) 



Since C is a symmetric positive definite operator, we can define C 1//2 and C x l 2 such that 
c i/2 c i/2 = c and c i/2 c -i/2 = j Then the definition of x* directly yields 

Res 2 {il>)= [ C- 1/2 (x* -Ce{<f> h )) : C 1/2 e{^)dx, 



and the Cauchy-Schwarz inequality finally leads to ( 1201) . 
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Theorem 4.5 (Reliability of the estimator) Under the assumption of corollary \4-S\ 
we have 

(er ) 2 < 4 C N 4 (- + ^(3 + 2 V3) + 4(/x + A)) V (t 2 + 4) llTfc " vIIt 
+4 (c + ^(3 + 2^+4^ + ^ + l) (^ + A)||C- 1 /2 (x *_ C£( ^ )) ||2 

+ max ((7 + 4V3); 2 + 2^ + ^(3 + 2^) + % + A)) (ji + X)c%j \\(j) h - R h (j) h \\ 2 H{ro 

+osc 2 (g). 

Proof: The theorem is a direct consequence of corollary I4.2[ lemma 14.31 and lemma 14.41 



Remark 4.6 In theorem 14. 5[ all constants are explicitly given. Indeed, even if cp and cr 
depend on the domain Q whereas k 2 and N depend on the used mesh, they can be evaluated 
or at least bounded, see [TU] and section below devoted to the numerical validations. 



Remark 4.7 The assumption t < a/ 3\c 2 f / Jl needed in corollary I4.2l is not restrictive since, 
in the Reissner-Mindlin model, t is expected to be a very small parameter, so that this 
property is naturally obtained. 

5 Efficiency of the estimator 

In order to prove the efficiency of the estimator, each part of it has now to be bounded 
by the error e r h ot up to a multiplicative constant. In the following, the notation a < b and 
a ~ b means the existence of positive constants c\ and C2, which are independent of the 
mesh size, of the plate thickness parameter t, of the quantities a and b under consideration 
and of the coefficients of the operators such that a < c 2 b and c\ b < a < c 2 b, respectively. 
The constants may in particular depend on the aspect ratio a of the mesh. 



Lemma 5.1 

( td ru n2 

r(rot,n) 



Rh - I)<Ph\\mrotn) < A- 2 * 4 ||7 - Tfclln + k - u h \{ 



Proof: Since 
we have 



+ |0-^| 2 + A^ 4 ||rot( 7 -7^ 
(R h - I)(j) h = A"V( 7 - 7 h ) - V(c^ - u h ) + (0 - <f> h ), 



\\R h - l)<f> h \\ < A t ||7 - 7 ft 1 1 + \uj - u h \i + ||0 - (fih\\, 
and with the Poincare-Friedrichs inequality, we get 



(Rft - /)^ir ;$ a-^ 4 ii 7 - 7/ill + k - + 
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Moreover, we have 

\\rot(<f> h - R h ^ h )\\ 2 < A- 2 t 4 ||rot( 7 - 7 ,)|| 2 + \<f>- <p h \l 
so that lemma 15.11 holds. ■ 

Lemma 5.2 There exists a relevant choice of x* such that 

Iic-^V - Cs{4>h))\\ 2 < hh - + 10, - <t>\l (21) 

Proof: First, there exists only one pair (0^,0*/*) £ Hq(Q) 2 x Q h solution of 

<M = -M) V^£# W, 
a((f)* h *,^ h ) = -(7&,Vfc) V4e6 fc . 
Then, by Theorem 3.9 of [30] and a relevant construction of a;*, for all T in 7/j we have 

iic-vv-c £ (0r))ii T < ii0*,-0riic. T . 

Because of the mesh regularity, we also get the global estimate 

iic- i /v-c £ (0r))ii<ii0^-0riic (22) 

Clearly 

C-V*( x * - Ce{<f> h )) = C- l ' 2 {x* - Ce(0D) + C x / 2 £(0r - 0ft)- 
By (122j) and the triangular inequality, we arrive at 

||C-VV " ^(0ft))ll < UC^V - Ce(C))ll + IIC -*h\\c 

<U* h -r h *\\c + K-M\c (23) 

Now, it remains to bound each of the two terms of the right-hand side of (|23|) . To begin 
with, let us consider iph £ 0^. Thanks to the definition of 0", we get 

a(0ft - 0h*> ^ft) = (7ft, *i>h ~ Rh^Ph) 

= {ihi iPh) ~ (7ft, Rh^h) 

= (j h - 7, if) h ) + a(4> h - 4>, tph) 

<(||7fc -7II-1+ |0ft -0|i)|^|i- 

By taking iph = 4>h — 0/1*) we obtain 

H0ft*-0ft||c< ||7ft- 7ll-i + 10ft- 0|i- (24) 
15 



Then, by the triangular inequality, we get 

\\4>l ~ M\c < \\<f>l ~ 0llc + ||0 - 4>h\\c, 
and by the definition of 0£, we have for all x/i G Hq(Q) 2 

so that 

HI ~ <P\\c < ||7-7fc||-i- 

We then obtain 

HI - <t>h\\c < \\l - lh\\-i + ||0 - M\c % 117 - 7h||-i + 10 - Mi- ( 25 ) 
Using (JMD and Qgft in flU, we get (J2TJ) . ■ 

Lemma 5.3 There exists a relevant choice of y* such that 

£ ^ + ^)ll7fc - Vlir £ ^Il7 - 7,f + ||7 - 7fc|| 2 -i + osc 2 (g), (26) 

T&T h 

where osc 2 (g) is an oscillating term. 

Proof: Because of lemma 3.1 of [15], we have for any T G Th the equivalence 

\\lh - y*\W ~ ^r 2 £ ll^ 7h - f*) ' ^1'^' 

where vt is the outward unit normal vector to T. Now we define y* as in [15], by noticing 
that (JSJ) implies that 

(7/» Vufc) = Vvft G W h , 

hence there exist fluxes qe G Vi(E), for all edges E such that 

/ ih-Vv h = / + / g-Tfft Vv h eVi(T), 

JT JT JdT 

where gx = gE^E^T-, v e being a fixed normal vector to E. According to the definition of 
the BDMi elements there then exists a unique y T G Vi(T) 2 such that 

y T -v E = gs V£ C T. 

Hence we define y* such that its restriction to each triangle T is equal to y T . According to 
its definition y* belongs to H div (Q) and moreover according to [32], we have 

divy* = -U h g. 
16 



Then by the use of theorem 6.2 from [T] and the mesh regularity we get 

Wlh - V*\\t ^ hr 2 ^ Wbfh ■ ve]e\\e + ^2 h T>\\divj h + g\\ T ', 

EcdT\dn T'Cut 

where [v] E denotes the jump of the quantity v through the edge E. Consequently 

T&T h TeT h EcdT\dn 

+ E E ^T'i^ 2 + h\)\\div^ h + gW^i 

T&T h T'Clot / 27 n 

< J2 h E (t 2 + h 2 E )\\[ lh ^ E ] E \\l 

EcdT\dn 

+ ^2h 2 T {t 2 + h 2 T )\\div^ h + gf T . 
TeT h 

Using the classical edge bubble functions as well as elementwise inverse estimates, it is 
proved in |11J . section 4.3 that : 

£ h E (t 2 + h 2 E )\\[ lh -v E ] E \\ 2 E < J2h 2 T (t 2 + h 2 T )\\g-Tl h g\\ 2 T 

Ee£(n)\an Ter h (28) 

+ Il7-7ft||-i + ^ 2 ||7-7ft|| 2 - 

Moreover, with the use of classical element bubble functions as well as elementwise inverse 
estimates, it is also proved in [11], section 4.1 that : 

J2 h 2 T {t 2 + h 2 T )\\div lh + U h g\\ 2 T < t 2 ||7-7/ l f + \\l-lhW-! 

T&T h 

+ Y, h T(h 2 T + t 2 )\\g-Tl h g\\ 2 T . 

T&T h 

Now, from ( 12 7p associated to the standard triangular inequality : 

\\diwy h + g\\ T < \\diwy h + U h g\\ T + \\g - Tl h g\\ T , 
the use of (USD and flU} leads to (J2BJ. 



(29) 



Theorem 5.4 (Efficiency of the estimator) There exists a relevant choice of x* and 
of y* such that 



*l|2 
T 



ACNk 2 J- + ^(3 + 2V3) + 4(/i + A)) J2 ^ + h T) hh-y 
\^ ^ J TeTh 

+ ( 2 C (- + ^(3 + 2^3) + 4(/i + A)) + 2] 2(/i + A) HC^'V - Ce{^ h ))f 
\ n fi J 

+ max ( (7 + 4^3); 2 + (- + ^(3 + 2^) + 4(/x + A))2(/x + ~X)c 2 R J \\<f) h - R h <f) h \\ 2 H 
\ fi fi J 

<{el ot ) 2 + osc 2 {g). 



(rot,Q) 
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Proof: The proof is a direct consequence of lemma 15.11 15.21 and 15.31 



6 Numerical validation 

Here we illustrate and validate our theoretical results by a simple computational example. 
Let Q be the unit square ]0, 1[ 2 . We consider the exact solution (u, <fi) in Q of the Reissner- 
Mindlin problem ([I])-© given by 



/ l-2x \ 
x 2 (l — x) 2 
l-2y 



exp 



x(l-x) 2/(1 — 2/) 



and 



with 



UJ 



\ y 2 (i-y) 2 J 

1 - (2/i + \)\~H 2 (a(x) + a(y))j exp 
, x Qz 4 - 12z 3 + 12z 2 - 



x{l - x) y(l - y) 
Qz + l 



a(z 



z\l - zY 

The corresponding scaled transverse loading function g is given by 

1 1 



g = [2ji + A) (c(x) + c(y) + 2 a(x) a(y)) exp 



with 



c(z) 



x(l-x) ?/(l-i/) 
120z 10 - 600z 9 + 1620^ 8 - 2880z 7 + 3504z 6 - 2952^ 5 



+ 



1708.2 4 - Q5Qz' 



z 8 (l - zf 
156z 2 -20^ + 1 



z 8 (l-z) s 

This analytical solution is extended by on dQ to obtain (u,<f>) G Hq(Q) x Hq(Q) 2 . Here 
we take t = 1/1024, A = 1, \i — 1 and A = 1. The meshes we use are uniform ones 
composed of n 2 squares, each of them being cut into 8 triangles as displayed on Figure [1] 
for n — 4. The refinement strategy is an uniform one so that the value of the mesh size h 
between two consecutive meshes is twice smaller. In order to validate the reliability of the 
estimator, we consider the "discrete error" given by 



-TOt 

"h,dis 



|2 + A -l t 2|| 7 _ Tfc ||2 + A-2t4|| rot ( 7 _ lh) ||2 + ||p fc7 _ lh \ 



where Phj stands for the piecewise Pi-discontinuous interpolation of 7 on the mesh l~h- 
This discrete error is defined by approximating the norm of 7 — 7^ arising in e£°* 

(see 0) by its discrete locally computable version defined by 



\\Phl-lh\ 



-i,h = SU P 

v h £W h 



\(Phl-lh,v h )\' 



(30) 
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Figure 1: Mesh level corresponding to n = 4 and h = y/2/8. 



The computation of \\Ph1 — Tft-I I— 1 ft, i s now an easy task and simply corresponds to the 
determination of the largest eigenvalue of a classical generalized finite dimensional eigen- 
value problem. In order to validate the reliability of the estimator according to theorem 
14.51 the error estimator is defined by 



(Vh) 



2 



4 C N 4 ( - + ^(3 + 2V3) + 4(/i + A)) £ (t 2 + h 2 T ) \\ lh -y*\\ 2 T 
+4 (( + ^(3 + 2V3) + 4(^ + A)^+l^ (fi + \)\\C^ 2 (x* -CsM)\\ 2 
+ max ((7 + 4^3); 2 + 2 ^~ + ^(3 + 2^3) + 4(/i + A)^ (ji + A)c^ - R^||^ (ro¥ 

and we plot on Figure [2] the evolution of the computed effectivity index f]h/e r ^ is versus 
h. Here, the values of x* as well as y* are respectively computed in the same manner as 
in [15] and [30], in order to obtain relevant choices as required by theorem 15.41 to ensure 
the efficiency of the estimator. Practically, some fluxes qe through the edges E of each 
triangle of the mesh are needed, and have to be computed by solving local linear problems. 
In fact, in our tests, these values are explicitely defined. For y*, we use Qe = {{jh • ^e}}, 
where {{7^ • ve}} denotes the averaged value on the triangles on each side of E of 7^ • u E 
evaluated at the middle of E. For x*, we use g% = Ylxetf(T){{C £ ( < l > h)}}( x ) l/ E^x- Here, 
{{Ce((f>h)}}(x) is the averaged value over the triangles surrouding the node x of the piece- 
wise constant function on each triangle Ce(<f)h), and X x stands for the classical Pi-Lagrange 
basis function associated with the node x. Moreoever, for the construction of x*, the Ar- 
gyris basis functions have to be used (see section 4 of [30J as well as [17J for the practical 
implementation) . 

From ([16]) we have £ = 1. The Poincare-Friedrichs constant ce is here equal to l/(y/2n) 

since Q is the unit square. Because of the kind of meshes used (see Figure []]), we have 

12 

N = 8 and k 2 = 1 H — -j=- (see annex 1772"]) . Finally, it can be proved [23] that on the unit 

V 27T 



square, cr < 2 J hence below we take this upper bound for cr (while it is conjec- 
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tured that cr 



27T 



7T 



tj, see [16]). As expected by the theory, it can be observed that the 



O estimator/ error 



Figure 2: rjh/e r ^ is versus h 



computed effectivity index is larger than one. Moreover, it converges towards a constant 
close to one when h goes towards zero, so that the proposed estimator is asymptotically 
exact. 
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7 Annexes 



7.1 Proof of (SJ) 

Let us consider v G C^°(Q) 2 . Two integrations by parts yield 



2 / |£(f)pdx = / \Vv\ z dx + / Vw(Vv) J (ix 
In Jn Jn 



/ |Vw| 2 <i£ + / \divv\ 2 dx 
Jn Jn 



> / \Vv\ z dx. 
Jn 

Hence by a density argument we obtain 

\\Vv\\ < V2 \\e(v)\\ Vd G Hq(£1) 2 . 



Then, we recall 
so that 



Ce{4>) = 2/i5(0) + ATr(e(0))X, 



C£(0)e(0)c?x 

n 

2/i / e(0)e(0)dx + A / Tr(e(0))Xe(0)dx 
Jn Jn 

2/i / |e(0)| 2 c/x + A / {Tre(<p)) 2 dx 
Jn Jn 

2 



>A*IIW. 



This proves ([3]). 



7.2 Evaluation of K2 for the triangulation of section 6 

With the definitions given above, let us consider z an affine function on lot, so that Jz = z 
on T. With v and v h defined in the proof of lemma 14.31 and the triangular inequality, we 
get 

\\V(v-v h )\\ T < \\V(v-z)\\ T + \\VJ(v-z)\\ T 

From (|T9l) . we get 

\\V(v-v h )\\ T < (l + Kl )\\V(v-z)\l T . 
Defining A = Vz and considering ip G Hq(Q) 2 , we have 

||V(u - v h )\\ T < (1 + ki)||Vi; - A\l T < (1 + «0 (UVu - V>IU + ||^ - AIU) . 
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Now, z is chosen such that 



A=-±- I <pdx. 

\Wt\ J uj t 

By Poincare inequality, there exists > 0, depending on the patch ut, such that 

U - MIt < C„ T h T ||V^|U V $ e # W- 

So, 

\\V(v-v h )\\ T < (l + Ki)||Vi;-^|| WT + (l + Ki)C r WT /ir||VV'|| WT 

< (l + K 1 )max{l;C^}(||Vt;-V'IU + MVV'IU). (31) 

V * ' 

= K 2 

Now, it remains to evaluate K\ as well as C^ T . Let rj z be the nodal basis associated to Wh- 
We have 

from what we deduce 

VJv = J2(v z -v)Vr] z ,Vv e Hl(tt). 

Let us define J\f T = J\f HT. We have 

||VJu||t = || ^ (v z - v)Vi1z\\t 

zeAfr 

- ^2 \\ Vz ~ v \\ t iiv^iit 

< Yl K-«IUI|v^|| T 

But 

||V?7 Z || T < p^ 1 , 

and from [TQl (5.12)], we get 

IK ~~ ^ c(w 2 ,2)||Vv|| Wz . 
With the triangulation involved, we have 

y/2hr 



and 



7T 



.„ T 3^2 h T „„ 
VJv T < — Vv 

7T p T 
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so that 

3v/2 hr 

K\ < . 

7T p T 

For the involved triangulation hx/pr — 2 and hence 



12 , N 

«i < (32) 

V2 7T 



3 

Since from [10], we have C WT = — , (ED and f[3"2"j) lead to 

12 

k 2 < 1 + 



V2 vr' 
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